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ABSTRACT 

We analyze the role of bars in the build-up of central mass concentrations in massive, disk 
galaxies. Our parent sample consists of 3757 face-on disk galaxies with redshifts between 
0.01 and 0.05, selected from the seventh Data Release of the Sloan Digital Sky Survey. 1555 
galaxies with bars are identified using position angle and ellipticity profiles of the /-band light. 
We compare the ratio of the specific star formation rate measured in the 1-3 kpc central region 
of the galaxy to that measured for the whole galaxy. Galaxies with strong bars have centrally 
enhanced star formation; the degree of enhancement depends primarily on the ellipticity of 
the bar, and not on the size of the bar or on the mass or structure of the host galaxy. The 
fraction of galaxies with strong bars is highest at stellar masses greater than 3 x 10 10 M G , stellar 
surface densities less than 3 x 10 8 M o and concentration indices less than 2.5. In this region of 
parameter space, galaxies with strong bars either have enhanced central star formation rates, 
or star formation that is suppressed compared to the mean. This suggests that bars may play 
a role in the eventual quenching of star formation in galaxies. Only 50% of galaxies with 
strongly concentrated star formation have strong bars, indicating that other processes such as 
galaxy interactions also induce central star-bursts. We also find that the ratio of the size of the 
bar to that of the disk depends mainly on the colour of the galaxy, suggesting that the growth 
and destruction of bars are regulated by gas accretion, as suggested by simulations. 
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1 INTRODUCTION 

It has long been conjectured that spiral galaxies evolve along 
the Hubble sequence from "late-type" disk-dominated systems 
to "early-type" galaxies with very massive bulges . This sce- 
nario is supported by the fact that the co-moving number den- 
sity of massive bulge-dominated galaxies increases from high red- 
shifts to the present day, while the number density of lower mass 
star-f orming disk galax i es evolves comparatively little with r ed- 
shift l lCoucn et al.lll998l ; lFasano et alj|200d : Ikovac et alj|20ich . It 
is believed that major merger events build classical bulges with 
light profiles that follow a r 1 ' 4 law dBarnes & Hernquistl 1 19921 : 
iToomrel fl977l) . while pseudo-bulges with light profiles that are 
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close to exponential are formed from the disk itself (see re- 
views bv lKormendv & Kennicuttll2004IWvse et alJI 19971) . Pseudo- 
bulges are common in disk-dominated galaxies at low redshif t 
jBalcells et alj2003l ; lLaurikainen et al.l2007l ; IWeinzirl et al.l2009r) . 



There are many mechanisms that may cause disks to become 
unstable and gas to flow towards the center of the galaxy, eventually 
forming a bulge. Tidal forces exerted by an interacting companion 
cause gas to lose angular momentum and to flow into the central 
region of the gala xy where it forms stars at an elevated rate (a so- 
called "starburst ")( lMihos & Hernquistl 19961 : Idi Matteo et al .120071 : 
ICox et al~ll2008l) . More than half of the galaxies with the highest 
central specific star formation rates (SS FR) in the loca l universe 
are interacting with a close companion dLi et al. Alterna- 
tively, asymmetric structures in the disk, including bars and spi- 
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ral arms exert torques that d rive gas inwards (see the review by 
iKormendv & Kennicutill2004) . 

The mechanisms by which bars influence the flow of gas in 
galaxies have been studied in considerable detail using simulations. 
Bars induce gravitational torques that drive the gas towards the 
leading end of the bar where it is compressed and shocked, and 
where it loses energy by radiative processes, causing it to flow to 
the center of the g alaxy l|Athanassoulalll992l: iKnapen et ail 1 2002 1: 
Pineret all 1 19951 : iRegan & Teubenl |2004 IZurita & Perezl bood: 



Sheth et all 120021). Galaxies with thick bars tend to form nuclear 
rings jRegan & Teubenl 120031 : IPineret al.lll995l: ISchinnerer et al.l 
120031) . So me of the gas will reach th e center and produce a disc- 
like bulge d Athanassouil 1 991 12005^ Bar s are robust if the galaxy 
remains undisturbed dAfhanassoulal Eo05bl) . 

Observational evidence that this basic picture is correct have 
come from highly resolved maps of the ionized, atomic and 
molecular gas in individual barred galaxies ( e.g. lOuillen et al .1 
Il995l : ISchinnerer et"aT1l2002t ISheth et al.ll200oh . which allow one 
to track the motion of the gas in detail. In addition, there 
have been studies of samples of barred and unbarred galaxies 
that have revealed h i gher gas concent rations djogee et al.l 120051 : 
ISakamotoet"aTir9 99; Shet h et all 20051) . higher cent ral star forma 
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1994H in the former. On average, these differen ces are most pro- 
nounced in galaxies with the strongest bars dHo & Filippenkol 
1 19971 : iMartiiJl 19951 : iMartin & Jean-Rene1ll994l) . Recent observa- 
tions suggest that most disk galaxies at high redshift are extremely 
gas-rich and hence dynamically unstable, and will likely form 
bulges on a short timescale (Genzel et al. 2006, 2008). 

It would be valuable to quantify the role that bars have played 
in the formation of the present-day bulge and pseudo-bulge popu- 
lation. It has been controversy about whether bars a re found more 
frequ e ntly in early-type or in late-type disk gala xies dBarazza et al .1 
120081: iMarinova et al'l [20091 : ISheth et al.ll2008l) . Enhanced central 
star formation rates are more likely to be found in early-type galax- 
ies with bars than in la te-type galaxies with bars dDevereuxl 19871 : 
lHo& Filippenkoli 19971) . These results are somewhat difficult to in- 
terpret because Hubble type is a complex function of both bulge- 
to-disk ratio and star formation rate. 

In semi-analytic models of galaxy formation, disks form when 
gas cools and settles at the center of a dark matter halo, while con- 
serving angular momentum. In recent years, simple prescriptions 
for bulge gro wth through disk instabilities have been included in 
these models dBenson et alj201ol : lDe Lucia et alj201 lMParrv et al .1 
120091) . but these have only been tested against observations in a 
rudimentary way. Ideally, these models should be able to account 
for the fraction of disk galaxies with pseudo-bulges as a function 
of both halo mass and stellar mass and the size of the disk. 

In this paper, we study the effect of bars in building the cen- 
tral mass concentrations in galaxies as a function of parameters 
such as stellar mass, stellar surface mass density and concentration. 
We ma ke use of imaging data from the Sloan Digital Sky Survey 
(SDSS. lYork et alj|2000h . The SDSS images have a pixel scale of 
0."396 and a mean PSF (point spread function) of l."4 (FWHM). 
In this study, we identify bars from the images of ~4000 disk galax- 
ies with redshifts in the range 0.01 to 0.05 drawn from the seventh 
data release (DR7) of the SDSS. We first investigate the bar strength 
needed to cause enhancement in the central star formation rate of 
the galaxy. We then divide our galaxies into a 4-dimensional pa- 
rameter space of stellar mass, concentration, stellar mass surface 



density and colour and quantify the fraction of galaxies that are 
undergoing bar-driven central mass growth at the present day. 

In Section 2, we describe our sample and the methods used to 
identify bars and to quantify their strength. In Section 3.1, we ex- 
amine how strongly star formation is concentrated towards the cen- 
ters of barred galaxies compared to control samples of non-barred 
galaxies matched in stellar mass, stellar mass surface density and 
colour. In this analysis, we parameterize bar strength using the el- 
lipticity of the bar, et ar - In Section 3.2 we examine the colours of 
bars, showing that the strongest bars have the bluest colours. In 
Section 3.3, we show how the fraction of barred galaxies varies ac- 
cording to the location of the galaxy in the 4-dimensional parame- 
ter space of stellar mass, stellar mass surface density, concentration 
and colour. In Section 3.4, we investigate which galaxy properties 
correlate most strongly with bar strength and size. Finally, in sec- 
tion 4, we summarize and discuss our results. 



2 DATA 



2.1 The Sample 



We select galaxies with M* > 10 10 M o in the redshift range 
0.01< z <0.05 from the MPA/JHU spectroscopic catalogue 
(http://www.mpa-garching.mpg.de/SDSS/DR7/), which is drawn 
from the seventh data release (DR7) of the SDSS dAbazaiian et al.l 
120091) . This yields a sample of 16573 galaxies. Our aim is to assess 
the role of bars in building bulges in present-day disk-dominated 
systems, so we also select systems with Rgo/Rso < 2.6, where R90 
and R50 are the radii of the circular apertures enclosing 90 percent 
and 50 percent of the total r-band light from the galaxy. We fur- 
ther limit the sample to galaxies with ellipticiy 1 - b/a < 0.25 
(or an inclination of less than 41.4 degrees), where a and b are the 
major and minor axes of an ellipsoidal fit in the r— band to each 
galaxy. We employ this criterion to select relatively face-on galax- 
ies, because it is difficult to i dentify bars in more inclined galax- 
ies dLaurikainen & Saloll2002h . These cuts yield a sample of 3890 
galaxies, which we will use to characterize bars. 

In this paper, we deribe global star formation rates by fitting 
the observed FUV, NUV, u, g, r, i and z band fluxes of the galaxy 
with a library of model SEDs dSaintonge et alj|201 ll : IWang et all 
1201 ll) . The model SEP library is generated u sing stellar population 
synthesis models dBruzual & Charlol 120031) and includes model 
"galaxies" spanning a range in metallicity, age, star formation his- 
tory and dust attenuation strength. The probability that the observed 
SED can be described by each model SED is calculated, and the 
SFR of the galaxy is then given as the probability-weighted SFR 
averaged over the whole model library. We also use the A v derived 
using this method to correct the global colours of galaxies for in- 
ternal reddening. 

We also make use of star formation rates estimated directly 
from the SDSS spectra, which are obtained through 3 arcsecond di- 
amteter fibres. These star formation rates, which are taken directly 
from the MPA/JHU catalogue, are obtained by fitting the measured 
emission line flux es to a set of photo-ioniz ation+ stellar population 
synthesis models dBrinchmann et al.l2004l) . 

Finally, we derive morphologica l parameters, such as asym- 
metry index for all the galaxies (see IWang et al.l [20 1 lL for more 
details). 
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2.2 Bar Identification 

We use position angle and ellipticity profiles to identify bars in 
galaxies. Similar methods have been e mployed in many past papers 
(e.g. iBarazza et all 120081: Ijogee et alj 12004 iKnapen et alj|2002l; 
Laine et ail 120021 ; IMenendez-Delmestre et alj |2007| ; ISheth et al .1 



20031 : IWhvteetaf. 2002; Wozniak etalll 19951) . 



We first use SExtractor to estimate galaxy sizes (Petrosian ra- 
dius, R 50 and R90), ellipticities and disk position angles. We mask 
neighboring sources in the vicinity of the main galaxy. We then run 
the IRAF.eJIipse task twice on the SDSS r— band images over a 
range of radii extending from twice the Petrosian radius down to a 
radius of 0.8 arcsec from the centre of the galaxy. 

The sampling radius in the fit increases exponentially with a 
step of 1.1, and the center of each ellipse is allowed to vary. Dur- 
ing the first run of the ellipse task, we set R90 as the starting ra- 
dius. During the second run, we set the starting radius to max(R 50 , 
1.5"). During both runs, the initial values of ellipticity and position 
angle are set equal to the values measured by SExtractor for the 
whole galaxy. The two profiles are then merged and bad points are 
discarded. Bad points are defined as points where the ellipse task 
failed to find a solution, the error on the ellipticity measurement is 
larger than 0.05, or the error on the position angle is larger than 5. 

The fit is said to have been successful if the merged profile 
covers the radius ranging from 1.4 arcsec to R 90 , and if it covers 
more than 60% of the original range in radii that served as input 
for the IRAF.ellipse task. We also visually examined all the fits 
and found that our procedure was successful in almost all cases. 
Failures occurred for 133 galaxies (less than ~ 3.5% of the sam- 
ple), mainly because the galaxies were strongly dis turbed or had 
extend ed l ow surface brightn ess disks. We note that lBarazza et alj 
d2008h and ljogee et alj | |2004|) obtained a similar failure rate in their 
own analysis of bars in galaxies in the Sloan Digital Sky Survey. 
The remaining 3757 galaxies are well fitted with ellipses through- 
out the entire disk, and our analysis of bars will be confined to 
these systems from now on. We do not deproject the images or the 
profiles before classifying a galaxy as barred or un-barred. Such 
de-projections are subject to considerable uncertainty, particularly 
for more bulge-dominated galaxies. Errors introduced by projec- 
tion effects will be discussed later. 

The criteria we use to classify a galax y as "barred" are sim- 
ilar to those adopted by Ijogee et alj d2004h . hereafter J04, with a 
few changes that are detailed below. J04 required that (1) within 
the bar, the ellipticity should rise to a maximum value greater than 
0.25, and that the position angle within the bar should remain con- 
stant to within 20 degrees; (2) at the end of the bar, the ellipticity 
should drop by more than 0. 1 and the position angle should change 
by more than 10 degrees. Figure [TJillustrates the ellipticity and po- 
sition angle profiles of a galaxy that is classified as barred according 
to these criteria. 

We imposed the following additional requirements after visu- 
ally inspecting all the fits : ( 1 ) the inner radius of the bar should be 
less than 0.5 times the outer radius of the bar, and ellipticity should 
increase monotonically as a function of radius from 0.33 R ouler to 
the end of the bar. This requirement improves the ability o f the al- 
orith m to identify bars in galaxies with prominent bulges iMartirJ 
1995b showed that the size of the bulge is on average about half 
the size of the bar). It also prevents complicated spiral structures 
and inner rings from being mistakenly identified as bars. (2) The 
outer radius of the bar should be less than 0.8 times the semi-major 
axis of the outermost fitted ellipse. This requirement ensures that 
edge-on, disk-dominated galaxies, where the inclination has been 



incorrectly measured, do not enter our sample of barred galaxies. 
(3) if criterion 1 adopted by J04 is satisfied, and the ellipticity drops 
by more than 0.25 at the end of the bar, but the position angle does 
not change by more than 10 degrees, we still classify the galaxy as 
barred. 

The final requirement ensures that we do not exclude galaxies 
which have a long bar surro unded by an aligned ring Examples 
of such systems are given in Gadotti & de Souza d2003l) . We also 
note that IMenendez-Delmestre et alj d2007l) estimated that the J04 
position angle criterion could introduce incompleteness at the 10% 
level. Criterion 3 resolves this problem for bars which have high 
ellipticity (e har >0.5), which are the focus of our study. 

As well as classifying galaxies into barred and non-barred sys- 
tems, we compute two quantitative parameters; 1) the bar ellipticity 
e&ar, defined as the value of the ellipticity at the end of the bar, and 
2) the relative bar size D bar ID disk , where D bar is the diameter of 
the bar and is is the diameter of the 25 mag/arcsec 2 isophote 
measured in the g-band. The median value of e bar is 0.47 and the 
median value of T) bar jD disk is 0.3 for the galaxies in this study. 

The most distant galaxies in our sample are at z ~ 0.05, where 
2.5 times the FWHM of the SDSS PSF (~ 3.5 arcsec) corresponds 
to a physical scale of 3.8 kpc. Note that a 3.8 kpc bar aligned along 
the minor axis of an inclined galaxy with 1 - b/a =0.25 has a de- 
projected physical size of 5 kpc, so our sample is complete for bars 
larger than 5 kpc. 

To demonstrate that our 5 kpc cut produces an unbiased sam- 
ple of barred galaxies, we plot the fraction of barred systems as a 
function of redshift and galaxy size in the left panel of Figure [2] If 
the barred galaxy sample suffers from incompleteness because we 
miss bars with smaller angular sizes, we might expect the bar frac- 
tion to drop at higher redshifts, particularly in small galaxies. We 
see, however, that bar fractions do not vary with redshift at fixed 
physical size. The right panel of Figure [2] shows the fraction of 
barred galaxies as a function of galaxy size and ellipticity. As can 
be seen, the bar fraction does not vary with ellipticity, so there is no 
clear inclination-dependent bias in our estimates of bar fraction, at 
least for values of less than 0.25. We will refer to the sample 
of galaxies with bars larger than 5 kpc as Bar Sample A. 

One might worry, however, that by adopting a fixed cut in the 
physical diameter of the bar, one would be able to identify a larger 
proportion of the bars present in large galaxies compared to smaller 
ones. Another possibility is to adopt a fixed cut on D bar jD iiA . We 
have examined the distribution of the diameters D15 of the galaxies 
in our sample. We find a lower limit of D 2 s of around 16 kpc, corre- 
sponding to the diameter of the most compact galaxies with stellar 
masses of 10 10 M G . We have therefore constructed a second sample 
of barred galaxies with D bar jD iiA > 0.3, which we will refer to as 
Bar Sample B. 

In total, we identify 1555 barred galaxies with D bm > 5 kpc 
( 41.4% of the original sample of 3757 galaxies), of which 827 
(2 2%) have e bar > 0.5. We n ote t hat a simila r fracti on was found 
by|M annova & Joged <2007h and llogee et all d2004h when a sim- 
ilar cut in bar size and bar ellipticity was applied. Bar Sample B 
contains 506 galaxies, and 363 of them have e bar > 0.5. In general, 
Bar Sample B is most robust for investigating how the fraction of 
bars depends on galaxy structural parameters, such as stellar sur- 
face mass density or concentration. Bar Sample A is most useful 
for investigating the fraction of central starbursts that were trig- 
gered by bar-induced inflows. 

Another issue that may potentially bias our results is that if the 
galaxy has a bulge, it may be more difficult to measure the shape 
of the bar by means of the ellipse-fitting technique. The bulge will 
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make the surface brightness contours of the galaxy rounder, caus- 
ing e bar to be under-estimated. As a result, some barred galaxies 
will then be missing from our sample. One way we can quantify 
this effect is through simulations, which are described in detail in 
Appendix A. The main conclusion of this study is that the iden- 
tification of "strong bars" with e bm > 0.5 is reasonably robust to 
the presence of a bulge. The interested reader is referred to the Ap- 
pendix for more details. In the rest of the paper, our conclusions 
will be based on the analysis of bars with e bar > 0.5. 

2.3 Control samples 

There have been many papers that have found that barred spiral 
galaxie s have higher central star formation rates than un-barred 
spirals JdeJong et aljl 984: Devereux 1987; Hawarde net aljl986t 
Ho & Filippenkd 1 19971: iMartin & Jean-Renel 1 1994 IPuxlev et all 
19881 : IZaritskv et aL.lll9940 . The enhancement is more pronounced 
for strongly barred galaxies than for weakly barred galaxie s 
dHo & FilipDenkolll997l ; [Martinlll995l ; lMartin & Jean-Renil994l) . 

It should be noted that the star formation rates in both barred 
and un-barred galaxies vary strongly as a function of their structural 
parameters, such as concentration and stellar mass surface density. 
If we wish to quantify the enhancement in central star formation 
rate accurately, we need to create control samples that are closely 
matched with the barred galaxy sample in these parameters. 

For each barred galaxy from Sample A, we search the parent 
sample for a matching galaxy (selected without regard to whether 
or not it has a bar) with stellar mass difference A log M, less than 
0.15, stellar surface mass density difference AlogyU„ less than 0.1, 
and g - i colour difference A(g - f) less than 0.1. 

2.4 Photometry of bars and pseudo-rings 

As described in Section 2.2, the outer end of a bar is defined as the 
point where the ellipticity profile reaches its maximum. From now 
on, we will refer to the region within the ellipse enclosing the end 
of the bar as t he "barred region" o f the galaxy. According to the 
simulations o f lAthanassouial l ll992l) . gas will flow from the end of 
the bar towards the center of the galaxy, following a curved path 
which is determined by the strength of the bar and the central mass 
density of the galaxy. The region of the galaxy that we define to be 
barred will cover the path of the inflowing gas until it condenses 
into a ring at the inner Lindblad resonance (see Figure 3). 

We also define an outer ellipse with ellipticity equal to that 
of the disk measured at its 25 mag arcsec" 1 isophote in the r- 
band and with major axis equal to l.2D bar . The "pseudo-ring" re- 
gion of the galaxy is the defined as the region between the el- 
lipse enclosing the barred region of the galaxy, and this outer el- 
lipse. The motivation for defining such a ring is that an ultra- 
harmonic re sonance ring in disk galax ie s is of t en found at the end 
of the bar jAthanassoula et al.l l2009ll bl l20ld; iRegan et al.l |2002| ; 
lRegan&Teubenll2003h . Gas will tend to accumulate in this ring 
and if the density reaches high enough values, stars will form. 
These outer resonance rings are usually aligned paralle l to the bar 
and may extend out to the corotation radius of the disk dButa et al.l 
2010; Sch warzlll98rl) . The "pseudo ring" region defined in this 
study is thus designed to cover the location of potential resonance 
rings in our galaxies. 

For each barred galaxy in our sample, we measure the SDSS 
g - i colour within the barred region and within the "pseudo ring" 
region. We also measure colours within the same regions in the cor- 
responding control galaxy, even though these systems often do not 



have bars see the right panel of figure 3). (Note that when we define 
the corresponding barred and pseudo-ring regions of the control 
galaxy, we scale the major and minor axes of both regions by the 
ratio D2s(control galaxy )/Z?25 (barred galaxy), so that the outer radii 
of the two systems match exactly^. We will compare the colours 
between the barred galaxies and control galaxies in Section 3.2. 

All the colours in the barred and 'pseudo ring" regions are 
corrected for dust attenuation using the average of the central at- 
tenuation and the global attenuation. The central attenuation is de- 
rived from the B aimer decre ment measured from the SDSS fibre 
spectrum. We have assumed a lCalzetti et al. I d200d) extinction curve 
and the line-to-continnum attenuation ratio from lwild et al ]l[20ll|). 
We apply a dust correction when both the Ho- and H/3 fluxes have 
S /N > 3. If the emission lines are too weak to measure the Balmer 
decrement, we adopt T v (star, fiber) from the measurements of the 
attenuation of the stellar continuum of the galaxy provided in the 
MP A/JHU catalog (tauv_cont ). This parameter is obtained by fit- 
ting HmzuaT&CMo3 d2003h population synthesis models to the 
stellar continuum; the reddening may then be estimated by deter- 
mining the extra "tilt" that must be applied to the models in order 
to fit the shape of the observed spectrum. The global attenuation is 
derived using SED-fitting techniques. 



3 RESULTS 

3.1 Establishing a link between strong bars and enhanced 
central star formation in galaxies 

As discussed in Section 2.3, we have constructed a control sam- 
ple with the same stellar masses, structural properties and global 
colours as the barred galaxies. By examining whether the central 
star formation rates are enhanced in the barred galaxy sample when 
compared to the controls, we can test whether there is a causal link 
between bars and central starbursts. By binning the barred galaxy 
sample and the control sample according to stellar mass, stellar sur- 
face density and morphological type, we ascertain to what extent 
the degree of central star formation enhancement is influenced by 
global galaxy properties. 

In this paper, we use the quantity sSFR/,/j/sSFR, ,, which we 
will denote C(SF), as a measure of how concentrated the star for- 
mation is in a galaxy: sSFR/vs is the specific SFR measured within 
the 3 arcsec SDSS fiber, and sSFR, D , is the global specific star for- 
mation rate of the galaxy. 

In Figure[4] we compare the distribution of log C(SF) for our 
barred galaxies (shown as black histograms on the plot) and con- 
trol galaxies (red dashed histograms). The two columns on the left 
show results for galaxies with weak bars (e bar < 0.5), while the 
columns on the right are for strong bars (e bar > 0.5). Each row 
in the figure shows results for galaxies split according to a set of 
global galaxy properties that include stellar mass, stellar surface 
mass density, concentration parameter, g - i colour, galaxy asym- 
metry index A, and ratio of the diameter of the bar to that of the 
disk (D b /D d ). In each case, the split is made at the median value 



1 We test two ways of measuring colour from a "barred region" (defined by 
the barred galaxy) in the control galaxy. In the first way, we do not change 
the orientation of a "barred region" and in the second way, we orientate 
the "barred region" so that its major axis points to the same direction as 
the major axis of the disk. These two types of colour measurements yield 
very similar results. In what follows, we present results using the first set of 
colour measurements. 
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Figure 1. An example of a galaxy identified as barred. In the left column, the top panel shows the image of the galaxy, the middle panel shows the fitted 
ellipses superimposed on the galaxy image, while the bottom panel shows the same image with the barred region highlighted. ^From top to bottom, the right 
column shows the surface brightness, ellipticity and position angle profiles of the galaxy. 
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Figure 2. The distribution of galaxies in our sample is plotted in the plane of physical size versus redshift (left), and physical size versus ellipticity (right). The 
coloured contours denote the fraction of bars as a function of position in this plane. Only Sample A barred galaxies are plotted. 
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of the parameter under investigation. The difference in the average 
log C(SF) and Kolmogorov-Smirnov probability for the barred and 
control galaxies for each sub-sample in each plot are sumarized in 
Table □ 

Our conclusion is very simple. For each of the parameters 
shown in the two left-hand panels (e bar < 0.5), there is no sig- 
nificant difference between the black and red histograms, i.e., there 
are no significant differences between the barred and the control 
samples. For each of the parameters shown in the two right-hand 
panels (e bar > 0.5), the differences ara very significant. This shows 
that only strong bars are capable of inducing enhanced central star 
formation in galaxies. In the right two panels, the offset between 
the black and red histograms is roughly the same for each parame- 
ter. This shows that the global galaxy properties considered here do 
not play a significant role in determining the degree of central star 
formation enhancement caused by bar-driven inflows. 

In the left panel of Figure [5] we show how the degree of 
central star formation enhancement depends on both e bar and on 
the concentration index of the galaxy (note that the concentration 
index is a simple measure of bulge-to-disk ratio). We plot con- 
tours of Ab- C logC(SF) in the plane of e bar versus R90/R50, where 
A fc _ c log C(SF) is the difference in logtsSFR^/sSFR,,,,] for the 
barred galaxies and the controls. We see that when e bar > 0.5, 
Ab-c logC(SF) is positive, indicating that star formation is more 
concentrated in the barred galaxies. When e bar < 0.5, star formation 
in barred galaxies and control galaxies are similarly concentrated. 

We note that there is an apparent deficit of galaxies with large 
concentrations and with high values of e bar in Figure [5] As dis- 
cussed in Appendix A, one possible reason for this is because e bar 
will be systematically underestimated in galaxies with larger bulge- 
to-disk-ratios. In Figure A2, we see that e bar is underestimated by 
about 0.1-0.15 when R 90 /R 50 > 2.5, which agrees well with the ob- 
served deficit. After taking this bias into account, we conclude that 
there is no clear evidence that strong bars are less efficient at chan- 
neling gas to the central regions of galaxies if a bulge component 
is already present. On the other hand, Figure [5] show that the de- 



gree of central star formation enhancement appears to be smaller 
in strongly barred galaxies with low central mass concentration 
(C < 1.9), in agreement with simulation results from dAthanassoulal 
Il992h and lShethetaHd2000l) . 

In the the middle panel of Figure [5] we investigate whether 
the degree of central star formation enhancement exhibits any de- 
pendence on the length of the bar relative to that of the disk. We 
plot A;,_ c log C(SF) in the plane of D bar jD dis k versus concentra- 
tion index. We see that A b _ c log C(SF) is somewhat higher when 
D bm ID disk > 0.3. We note however, that this effect is mainly due to 
the fact that the size of the bar is correlated with its ellipticity. This 
is clearly demonstrated in the right panel of the figure, which shows 
contours of A b _ c log C(SF) in the D bar /D disk versus e bar plane. The 
main dependence of C(SF) is on e bar and not on D^/D^. 

In summary, we find that the central star formation enhance- 
ment depends mainly on the ellipticity of the bar, and not on the 
size of the bar or on the mass or structure of the host galaxy. 



3.2 Bar and pseudo-ring colours 

In this section, we compare the colours of the barred and pseudo 
ring regions in the barred and control galaxy samples. As shown in 
the previous section, the degree of central star formation depends 
prinarily on ellipticity, so we simply split the sample at e bar = 0.5. 
The left panels of Figure [6] show that the colours in the barred and 
pseudo-ring regions do not differ significantly from those measured 
in the control samples when e bar < 0.5. Much stronger differences 
are found when e bal . > 0.5. These results imply that the enhanced 
central star formation extends well beyond the central region of the 
galaxy when the bar is strong. Enhanced star formation is found in 
regions where gas is likely to flow inwards, as well as in the outer 
regions of the galaxy, where gas is compressed by the resonance. 

Finally, in the left panel of Figure [7] we examine how the av- 
erage ellipticity of the bar varies as a function of location in the 
plane of bar colour versus global galaxy colour. As can be seen, if 
the barred region of the galaxy is blue, then a strong bar is gen- 
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Figure 4. Histograms of the ratio of specific SFR evaluated within the SDSS fiber to specific SFR for the galaxy as a whole (log C(SF)) are shown for 
the barred galaxy sample (open histograms), and for the control sample matched in stellar mass, stellar mass surface density and g - i colour (red, hatched 
histograms). The plots on the left show results for galaxies with weak bars {e^ ar < 0.5), while the plots on the right are for galaxies with strong bars (£W > 0.5). 
The samples are further divided by stellar mass (log M», M» is in units of M Q ), stellar mass surface density (log//,, fi, is in units of Mo/kpc?), concentration 
(R90/R50), g - i colour, asymmetry index (A) and relative bar size (Pbar I Ddisld at the median value of each parameter, as denoted in each plot. Only Bar Sample 
A galaxies are plotted. 
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Figure 5. The distribution of C(SF) of barred galaxies minus C(SF) of control sample galaxies is plotted in the plane of bar ellipticity parameter ei, ar versus 
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and Di, ar of their corresponding barred galaxies when making the plots. Only Bar Sample A galaxies are considered. 
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Table 1. The difference of the C(SF) between the barred and control galaxies (see Figure EJ. Aj,-c log C(S F) is the average value of logC(SF,barred)- 
logC(SF,control). KS prob is the Kolmogorov-Smirnov probability that the two distributions (C(SF) for the barred and control galaxies) are drawn from 
an identical parent population. The columns on the left show results for galaxies with weak bars (eb ar < 0.5), while the columns on the right are for galaxies 
with strong bars (ei, ar > 0.5). The samples are further divided by stellar mass (log M„, M„ is in units of M Q ), stellar mass surface density (\og fi„, fi, is in units 
of Mo/kpc 2 ), concentration (Rgo/RjQ), g — i colour, asymmetry index (A) and relative bar size (Db ar l F>disk) at the median value of each parameter, as denoted 
in the second row of the table. Only Bar Sample A galaxies are calculated. 



erally located within it (note that the "noise" in the plot probably 
arises because the barred region only traces the gas inflow region 
in an approximate way). On the other hand, there is no obvious ef- 
fect seen as a function of the globed colour of the galaxy. In the 
right panel, we examine how Oior/Afcfc varies as a function of the 
same two parameters. In contrast to the results obtained for eb ar , the 
size of the bar (scaled to that of the disk) does not correlate with 
the colour of the bar, but scales strongly as a function of the global 
colour of the galaxy. The reddest galaxies host the longest bars. We 
will discuss possible explanations for this result in Section 4. 



3.3 Fraction of galaxies with strong bars 

In this section, we examine how the fraction of disk galaxies with 
strong bars (i.e. galaxies having ei, ar > 0.5) varies as a function 
of stellar mass, stellar mass surface density, and concentration. We 
show results for bar samples A and B. Sample A, which includes 
all galaxies with bars larger than a physical size of 5 kpc, is use- 
ful for deriving a stronger lower limit to the fraction of galaxies 
with strong bars. Sample B, which include only those galaxies with 
Dbarl Djjsk > 0.3, should yield unbiased trends in strong bar frac- 
tion as a function of parameters such as stellar surface mass density, 
which scale with the size of the galaxy. 

In the top row of Figure[8] we plot contours of strong bar frac- 
tion from Bar Sample A in the plane of stellar mass versus stellar 



surface mass density (left), and stellar mass versus concentration 
index (right). The fraction of galaxies with strong bars exceeds 
0.5 for galaxies with logM,/M G > 10.6 with low stellar mass 
surface densities (log/j„/M G fcpc~ 2 < 8.5) and low concentrations 
(C < 2.5). For lower mass galaxies, the strong bar fraction is a fac- 
tor of two smaller, but also peaks at low stellar surface densities 
and concentrations. 

In the bottom row of Figure [8] we show results for sample 
B. The strong bar fractions are lower and the plots are noisier, but 
the same qualitative trends are apparent. We conclude, therefore, 
that the bar-driven mode of bulge formation appears to be most 
ubiquitous in massive, disk-dominated galaxies at the present day. 

Finally, in Figure [9] we plot contours of strong bar fraction in 
the plane of C(SF) versus stellar mass, stellar surface mass den- 
sity, and concentration parameter. This plot should tell us whether 
strong bars are not only a sufficient condition for enhanced central 
star formation in galaxies, but whether they are also necessary con- 
dition. Note that in this plot, A(log C(SF)) is defined as log C(SF) 
minus the average value for all galaxies with similar stellar mass 
mass surface density and concentration from the "parent sample": 
positive values indicate that star formation is more concentrated 
than average, negative values that it is less concentrated. 

We see from the top row of the figure that the strong bar frac- 
tion peaks in galaxies with the most concentrated star formation 
and with stellar masses greater than ~ 3 x 10 10 M o , stellar surface 
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Figure 6. The top panels show a comparison of the distribution of g ~ i colours measured for the " barred regions" of the galaxies classified as barred (black 
histograms) with the distribution of g - colours measured for the same regions in the sample of control galaxies (red dashed histograms). Results are shown 
for weak bars (left), and for strong bars (right). The bottom panels show a comparison of g - i colours measured for the "pseudo-ring" regions of the barred 
and control galaxy samples. Only Bar Sample A galaxies are plotted. 



densities less than ~ 3x 10 8 M o kpcT 2 , and concentration parameters 
less than ~ 2.4. The peak strong bar fraction is around 0.5. This im- 
plies that although bar-driven inflows constitute a major channel for 
inducing enhanced central star formation in disk-dominate d galax- 
ies, th e y are not the only pro cess at work. As shown by iLi et al.l 
uReichard et al. galaxies with close companions and 

lopsided galaxies also have enhanced star formation rates. We note 
that in the highest surface density galaxies, bars appear to play no 
role at all in inducing enhanced central star formation. 

Very interestingly, there appears to be a secondary peak in the 
fraction of strong bars in galaxies with the same stellar masses, 
stellar surface densities and concentrations, but with central star 
formation rates that are significantly below average. In other words, 
there also appears to be a population of quenched galaxies with 



strong bars. This is consistent with the finding in lSheth et alJfeOOsT) 
that some barred spirals have nuclear regions that are deficient in 
gas. 

This result raises the question as to whether central starbursts 
produced by bar-driven inflows play any role in shutting down star 
formation in some galaxies, either by using up the available gas 
more quickly, or by generating winds/outflows that expel the gas 
from the central regions of these systems. We find no enhancement 
of optical or radio AGN fraction in this population, which suggests 
that starburst-driven processes may be at work. Another possible 
explanation is that the gas has simply been consumed and that the 
galaxies are in temporarily quiescent state. This will be the subject 
of a future investigation. 

Finally, we note that for the low density, low concentration 
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Figure 7. The barred galaxy sample is plotted in the 2D plane of bar colour versus global galaxy colours. In the left panel, coloured contours denote the 
average e bar as a function of position in the plane. In the right panel, the contours indicate the average value of D bar ID c n s k- Only Bar Sample A galaxies are 
plotted. 



galaxy population < 3 x 10 M e kpc , C < 2.5), the strong bar 
fraction is smallest when there is no enhancement in central star 
formation rate. This is consistent with the idea that in the absence 
of any dynamical perturbations, the star formation rates in present- 
day disk galaxies are regulated by accretion from the surrounding 
halo (e.g. Kauffmann et al 1993). Bar-driven inflows then play a key 
role in regulating the rate at which the accreted gas in consumed 
into stars. 



3.4 The relation between bar properties and global galaxy 
properties 

In Figure |10|and Figure [TT] we examine how e bar and D bar jD iisk 
depend on global galaxy properties such as stellar mass, concentra- 
tion index and g - i colour. 

The left panels of Figure [TO] show that more massive and less 
concentrated galaxes tend to have stronger bars. This is in agree- 
ment with results discussed previously, showing that the fraction of 
galaxies with strong bars peaks for massive, disk-dominated galax- 
ies (see Appendix for a discussion about bulge contamination). The 
right panel shows that at a fixed concentration, there does not ap- 
pear to be any direct connection between e bar and global g—i colour. 
Once again, results for Bar Sample A and B are qualitatively very 
similar. 

In Figure [TT] we again see that £W/A;,- S ,t is strongly corre- 
lated with global g - i colour. At a fixed stellar mass or concen- 
tration, redder galaxies have higher bar-to-disk size ratios. In some 
extreme cases, the whole galaxy consists of a long bar and a sur- 
rounding ring. These galaxies are almost all red. Figure 12 shows a 
compendium of SDSS images of such galaxies. It is interesting to 
note that at a fixed colour, D^/D^ does not vary as a function of 
concentration or stellar mass. 

In summary, e bar is correlated with both stellar mass and struc- 
tural parameters such as concentration, but not with global colour. 
Di,„ r /D djsk is very strongly correlated with global colour, and only 
weakly correlated with mass and structural parameters. 



4 SUMMARY AND DISCUSSION 

We have identified 1555 galaxies with bars longer than 5 kpc in a 
sample of 3757 face-on disk galaxies with stellar masses greater 
than 1O 1O M and redshifts in the range 0.01 < z < 0.05 drawn from 
the seventh data release of the Sloan Digital Digital Sky Survey. We 
have measured the ellipticity, length and g - i colour of the bars and 
have explored the relationship between global galaxy properties, 
bar properties, and enhanced central star formation in galaxies. 
Our main conclusions may be summarized as follows: 

(i) Only strong bars with ellipticity greater than 0.5 result in en- 
hanced central star formation in galaxies. This result holds for all 
galaxies, regardless of their stellar mass, stellar surface mass den- 
sity or bulge-to-disk ratio. (Section 3.1, Figure [4) 

(ii) Strong bars tend to have blue colours, and are surrounded by 
a blue outer ring. This indicates that the enhanced star formation is 
occurring on the scale of the bar itself. (Section 3.2, Figure[6} 

(iii) The incidence of strong bars is highest for massive galax- 
ies with low surface mass densities and concentrations, reaching 
values as high as 60% for galaxies with stellar masses greater than 
3 X 10 10 M o , stellar surface densities less than 3 x 10 8 M G kpc~ 2 md 
concentration indices less than ~ 2.5. (Section 3.3, Figure[8} 

(iv) Within this population, the incidence of strong bars appears 
to be bimodal. The fraction of bars is highest for galaxies with cen- 
tral star formation rates that are factor of 2-3 higher than the aver- 
age. There is a secondary peak in strong bar fraction for galaxies 
with central star formation rates that are more than a factor of 10 
lower than average. (Section 3.3, Figure[9) 

(v) There is no correlation between the ellipticity (or strength) 
of the bar and the global colour of the galaxy. There is a strong 
correlation between the size of the bar and the global colour of the 
galaxy. 



4.1 Relation of this work to other recent studies of bars in 
large galaxy samples 

There have been a number of studies focusing on how the frac- 
tion and properties of bars varies as a function of Hubble type, or 
parameters such as stellar mass. Because most galaxy properties 
are strongly correlated, it is difficult to identify which is the pri- 
mary parameter in many of these correlations. Our approach of 
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Figure 8. The distribution of galaxies from Bar Sample A (top row) and Bar Sample B (bottom row) is plotted in the 2-dimensional plane of stellar mass 
versus mass surface density and stellar mass versus concentration. The coloured contours denote the strong bar fraction as a function of position in the plane. 
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Figure 9. The distribution of galaxies from Bar Sample A is plotted in the 2-simensional plane of of A C(SF) versus stellar mass, A C(SF) versus mass surface 
density and A C(SF) versus concentration. The coloured contours denote the strong bar fraction as a function of position in the plane. 



analyzing the properties of bars in two-dimensional projections of 
different galaxy parameters helps us differentiate primary from sec- 
ondary correlations. 

Marin ova et al (M09) found an increasing bar fraction 

as a functio n of galaxy luminosit y and stellar mass, consistent with 
our results. iBarazza et al. 1 d2008l) (B08) found the opposite trend. 
However, our trends in bar fraction as a function of concentration 
index are consistent with B08. but at odds with M09. The main 
reasons for these discrep ancies are the very diff erent sample defi- 
nitions (see discussion in lNair & Abrahamll2010h . We note that the 
analysis in this paper has been carried out for a volume-limited 
sample of galaxies with stellar masses greater than 10 10 M G , so the 



statistics presented in this paper should not suffer from any sample 
selection biases. 



lLaurikainen et al.l d2007h found th a t ei, ar does not vary as a 
functio n of Hubble type. IWhvte et alj d2002t) and lAguerri et al.l 
d2009h found a tende ncy for later-type galaxies to have higher e 6ar , 
while iMartinl i ll 9951) found a tendency for early-type galaxies to 
have higher e bar . Our result show a possible explanation for these 
discrepancy from different studies: early-type galaxies are both 
more massive and more concentrated than late- type galaxies, while 
we show that the trend of ei, ar varying with stellar mass is contrary 
to the trend with concentration. We have chosen to bypass Hubble 
type altogether, focusing instead on conecntration index and stellar 
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Figure 10. Galaxies from Bar Sample A (top row) and Bar Sample B (bottom row) are plotted in the 2-dimensional plane of concentration versus stellar mass 
and concentration versus g - i colour. The coloured contours denote the average value of £(,„,. as a function of position in the plane. 
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Figure 11. Galaxies from Bar Sample A are plotted in the 2-dimensional plane of stellar mass and concentration versus g — i colour. The coloured contours 
denote the average value of Djor/A&i as a function of position in the plane. 



surface densities as parameters that have a simpler physical inter- 
pretation. 

There have been many attempts to link the length of the bar 
length with the global structure of galaxies. Longer bars were found 
to correlate with more prominent bulges in early - type galaxies 
( Aguerri et al. 2009; Athanassoula &Martinet 1980; Athanassoula 
120031 : lErwiijl2005l : lMartinlll995h . iHovle et al.l l l201ll) was the first 
to find that galaxies hosting longer bars also tend to be redder. Our 



results show that D bar /D djsk has the tightest correlation with colour 
and at a fixed colour Dt ar ID c n s k does not vary with stellar mass or 
concentration. 

There have also been numerous studies of the relation between 
bars, gas concentration, central star formation rate and metallicity 
gradients in galaxie s, supporting the picture that bars induce gas 
inflows in galaxies dde Jong et al.ll 19841; Ijogee et alj|2005t iMarthi 



ll995l ; lMartin & Jean-Rene1l994l;ISakamoto et al.ll999l ; ISheth e~ 
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Figure 12. A gallery of images of galaxies with long bars. 



120051 : IZaritskvetallll994) . Recently, lEllison et al.l d201 lh used a 
large, visually classified sample to show that barred galaxies with 
stellar M, > 10 10 M G have centrally enhanced star formation rates. 
This is consistent with our results. No dependence o f central star 
forma tion enhancement on bar ellipticity was found in lEllison et al.l 
J20T lj). One possible reason for this discrepancy is that their control 
sample was only matched in stellar mass, while our control sample 
was matched in stellar mass, surface mass density and colour. Be- 
cause the central star formation rates of galaxie s depend on a l l these 
parameters (and on colour in particular), the lEllison et al. I d201lh 
might not have had the sensitivity to isolate the ellipticity depen- 
dence. 

Finally, we note that our results are inconsiste nt with some 
past findings dDevereuxll987l : lHo & Filippenkoll997h that late type 
galaxies do not have significantly enhanced central star formation . 
One possible reason for this discrepancy is that these studies did not 
distinguish strongly barred galaxies from weakly barred galaxies. 
As we have seen, the fraction of galaxies with strong bars peaks 
at higher stellar masses, so the late-type samples investigated in 
these older studies may have only contained galaxies with weak 
bars. We also caution that bar classifications in bulge-dominated 
galaxies are subject to strong biases. Two dimensional image de- 
composition techniques that allow the bulge component to be prop- 
erly subtracted are one obv ious way forward (for example, iGadottil 
120091 : IWeinzirl et ai]|2009l) . but this lies well beyond the scope of 
this study. 

4.2 Comments on the results in this paper 

4.2.1 Only bars with e bar > 0.5 induce enhanced central star 
formation in galaxies 

We found that only bars with e bm > 0.5 were able to induce en- 
hanced levels of star formation in the central and bar-dominated 
regions of the galaxies in our sample. In addition, star formation 
in the circular region surrounding the bar was also found to be en- 
hanced. 

These results are qualitatively consistent with simulations that 
show that gas inflows are stronger in galaxies with strong bars, be- 
cause such bars induce stronger torques that pull the gas off their 
initially circular orbits, causing it to flow towards the center of the 
galaxy. Some of the in flowing gas will stop outside the corotatio n 
radius and form a ring dAthanassouialll992tlFriedli & Benzlll993l) . 

Theoretically, the torques induced by the bar are proportional 
to the mass and to the elongation of th e bar, but inversely propor - 
tional to the central axisymmetric mass dCombes & Sandersll98ll) . 



Observationally, these three factors correspond to the mass of the 
bar relative to that of the disk, the bar ellipticity, and the bulge to 
disk ratio of the galaxy. In order to determine the mass of the bar, 
we would need to subtract a model for the bulge and the disk of 
the galaxy, which is beyond the scope of this work. In this study, 
ellipticity has been used as our measure of bar strengt h, because it 



is the more easily accessible observational parameter ( Block et al. 
l200ll ; lLaurikainen et al.l2002l ; lMartinez-Valpuesta et al.ll2006t) .~ 



As the bulg e-to-disk ratio increases, the effect of the ba r will 
become weaker dLaurikainen et alj|2007l IWhvte et al1l2002l) . Pa- 
rameterizing bar strength in terms of e bar for all galaxies, inde- 
pendent of their bulge-to-disk ratio, does not take this effect into 
account. However, as we have discussed, e bar will be underesti- 
mated by our ellipse-fitting method in bulge-dominated galaxies. 
These two effects may somehow conspire to make e bar an excel- 
lent predictor of central star formation enhancement, but we admit 
that further simulation studies are required to make a more detailed 
interpretation of our results 



4.2.2 The fraction of strong bars in galaxies with central 
starbursts is peaked at value of 0.5 — 0.6 

Tidal interactions are another obvious mechan ism for inducing 
enhanced central star format i on in galaxies dCox et al. 1 120081 - 
iBournaud et al.l2007h . [Li et al.l J2008h found that >~ 40% of SDSS 
galaxies with the highest central specific star formation rates esti- 
mate d have a close companion within a projected radius of 100 
kpc). iLi et alj d2009h showe d that barred gal axies do not have an 
excess of close companions iMarinova et al .1 (simi lar results were 
also found by also found by |2009|);lBarazza et al.l (similar results 
were also found by also found by 2009). In our study, we have ex- 
cluded interacting and merging galaxies from the "parent sample" 
in order to ensure that the ellipse fitting procedure delivers robust 
results. 

In summary, our main conclusion seems to be that bars and 
galaxy-galaxy interactions are independent phenomena that to- 
gether could explain most of the central starbursts in galaxies in 
the local Universe. 



4.2.3 The size of the bar is strongly correlated with the global 
colour of the host galaxy 

Simulation results indicate that once a bar is formed, it will rapidly 
slow down and grow longer as angular momentum is transfered 
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to the bulge and to the halo. On the other hand, if gas con- 
tinues to accrete onto the galaxy, the gas may transfer angu- 
lar momentum to the bar, thereby re-accelerating it. This sup- 
presses that growth of the bar and can even cause it to be- 



come shorter ( Athanassoula 20021 20031 




Berentzen et all 20071; 


Curir et al.l 120071 Debattista & Sellwood 




20001: iFriedli & Benz! 


19931; Villa- Vareas et al.l 20ld: Weinberg! 


1985). Given that there 



is a tight correlation betw een colour and HI gas fraction in galax- 
ies dCatinella et al]|2010h . our results suggest that the observation 
effect may come about because the size growth of the bar is being 
regulated by gas accretion in galaxies. Red galaxies have not had a 
gas accretion event for many Gyr, so their bars have been able to 
grow. The fact that the colour of the bar and the size of the bar do 
not correlate, may imply that size of the bar depends on the gas ac- 
cretion rate integrated over the past history of the galaxy, whereas 
the colour of the bar depends on the current gas inflow rate. 



4.2.4 The incidence of bars peaks in disk-dominated galaxies 
with high masses and low stellar surface densities 

Galaxy formation models in a ACDM cosmology predict that disk 
instabilities represents the dominant contribution to the formation 
of bulges in massive galaxies with M, between 10 10 and 10" M Q 
dDe Lucia et al.|[2oTTh . In these studies, bars are assumed to form 
when disks satisfy the standard Toomre Q-par ameter criterion, 
which has been validated by N-body simulations dEfstafhiou et al.l 
[19821) : 



V m 



(GM imk /R d yP- 



<e~ 0.5-1. 



(1) 



We find that fjbar) peaks when M, is high, in agreement 
with the |Pe Lucia et alj 120111) results. However, at fixed stel- 
lar mass, f(bar) increases at lower stellar lower surface densities, 
When p > 10 8 5 M o kpc" 1 or R 90 /Rso > 2.6, f(bar) quickly drops 
to values near 0. This is not easily understood in the context of 
the formula given abo ve. Athanassoula (2008) has criticized the 
lEfstathiouetaild 19821) criterion, pointing out that it is not of gen- 
eral applicability and should not be used in semi-analytic models. 
In particular, stellar disks containing gas would be expected to be 
significantly more unstable than predicted. 

Studies have shown that p., ~ 10 85 M o kpc~ 2 is a spe- 
cial transition point for galaxies in the local Universe. At sur- 
face densities lower than this value, galaxies contain significant 
amounts of gas, but at higher surface densities, the atomic and 
molecular gas fraction in galaxies drop sudd enly in many galaxies 
dCatinella et ai]|20ld: ISaintonge et al]|201 lh . A similar transition 
is seen at R90/R50 ~ 2.6. In contrast, the correlation between gas 
content and stellar mass is not as pronounced. It is thus reasonable 
to suppose that the fraction of bars may peak at stellar surface den- 
sities below 10 85 M Q kpc~ 2 in the local Universe, because this cur- 
rently represents some special instability point for galactic disks. 
What we are not able to say from these observations, is what role 
bars have played over the whole history of the Universe in shap- 
ing this instability point. In any case, our results do suggest that 
gas accretion, disk instabilities, structural changes, and quenching 
of gas and star formation may happen coherently in galaxies. In 
the present-day Universe, bars are clearly an important regulatory 
mechanism in the "transition regime" between the blue and the red 
populations. 
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APPENDIX A: TESTS OF OUR BAR IDENTIFICATION 
PROCEDURE 

In Sect 2.2, we discussed that the existence of prominent bulges 
will make the elliptical isophotes of galactic bars rounder, causing 
us to miss bars in bulge-dominated galaxies and under-estimate the 
bar ellipticity. Here we quantify the affect of bulge contamination 
on our estimates of ellipticity and bar fraction by simulating the 
effect of the bulge of our measurements . 

First, we select the galaxies with Rgo/Rso < 2. 2 from the orig - 
inal barred sample. According to the formula from iGadottil (120091) . 
R90/R50 < 2.2 corresponds to a bulge mass fraction of less than 
14%. We use the images of these galaxies as "inputs", and we 
then add idealized bulge components onto them as described be- 
low. ^From now on, the bar parameters related to the original im- 
ages will be referred to as intrinsic and the parameters related to 
the simulated images will be referred to as the output. 

A pseudo bulge is more compact than the classical bulge 
dGadottil l2009h . and is thus is more likely to distort the inner 
isophote of galaxies, so to be conservative, we only simulate 
pseudo-bulges in this test. We set the stellar mass of the pseudo 
bulges to range from 0.03 to 1.8 times the stellar mass of the orig- 
inal galaxies. This range is chosen to produce the R ga /R 50 range of 
the galaxies in our "parent sample". 

^From Figure 13 of IGadottil d2009l) . we estimate log r e ~ 
0.2x log M, - 2.05, where r e and M, are the effective radius 
and stellar mass of the pseudo-bulge. We describe the surface 
brightness of a pseudo-bulge by a sersic profile with index 1.4, 
I(r) = Iq x exp(-2.47 x (r/r e ) 1/1A ), where I(r) is the surface 
brightness at radius r, and I p is the central surface brightness 
dKormendv & Kennicuttll2004l) . We assume that the /-band flux is 
linearly proportional to the stellar mass dBell et al.l2003l) . Integrat- 
ing the surface brightness profile to get the total flux, I is estimated 
as/ = fluxl(\.\l xr]). 

We produce 2000 simulated galaxy images by adding the sim- 
ulated pseudo-bulge components to the real /-band images of the 
disk-dominated galaxies. The distribution of concentration indices 
(Rgo/Rso) of the simulated galaxies is shown in Figure IaTI We can 
see that the simulated Rgo/Rso ranges from 1.8 to ~3, and follows 
th e same relation with bulge-to-total mass ratio of the galaxy shown 
in lGadottil J2009h . 

We then use the ellipse-fitting method described in Sect. 2.2 to 
identify and measure bars from the simulated images. We refer to 
the fraction of bars identified from the simulated images as the out- 
put bar fraction {F mnpul (bar)), and the correspondingly measured 
bar ellipticity as e bal%outpl „. 

In the left panel of Figure lA2l we show how the intrinsic ellip- 
ticity of the bar {ebarjntrinsic) is under-estimated as the concentration 
increases. We see that when Rgo/Rso =1.8, 2.4 and 2.8, e bar j„ tr j„ sic is 
under-estimated by 0.065, 0.15 and 0.18. When Rgo/Rso > 2.7, the 
error in flattens at a value of ~0. 18. 

In the middle panel of Figure IA2I we show how the e bart0Ulput 
from simulation and the e barobs measured from real galaxies vary as 
a function of concentration R 90 /Rso- The relation between e baroutput 
and Rgo/Rso (black line) shows the pure bulge contamination effect. 
We can see that the bulge contamination could be responsible for 
all the decreasing trend of e bar as a function of R 90 /Rso- Thus it 
is unclear whether there is a intrinsic trend that when Rgo/Rso is 
lower, e barintrinsic is higher. 

In the right panel of Figure IA2I we show the bar verification 
rate (F(bar)) from the simulated images in the plain of output con- 
centration Rgo/Rso versus input e bar j mri „ sic . We can see that at a 



fixed Rgo/Rso, F{bar) is strongly correlated with input e bar j ntr i„ S i C . 
At Rgo/Rso ~ 2.6, when e bar j mrimii: ~ 0.3, only ~ 20% of the 
barred galaxies are identified. This is expected because we require 
e ba r,output > 0.25 for a galaxy to be classified as barred. When 
e ba r,mtrinsic ~ 0.5, ~ 50% of the barred galaxies are identified. When 
e bar,mtinsic > 0.6, F(bar) flattens at a value of about 0.7. Bulge con- 
tamination thus results in loss of bars with low values of e bm . 

In this study, we focus on galaxies which have e barj , bs ~ 0.5. 
This corresponds to an input e bar in , rinsic that is a weakly rising func- 
tion of Rgo/Rso, as shown by the grey curve in the right panel of 
Figure lA2l . We see that most of the region above the e barMbs ~ 0.5 
curve corresponds to F(bar) values greater than 70%. In the region 
enclosed by the curve and e bar j,„ ri „ sic ~ 0.5, Fibaf) can vary from 
values of around 70% for the less concentrated galaxies to values 
of around 50% for the more concentrated galaxies. We conclude, 
therefore, that we are likely correctly identifying between 2/3 and 
3/4 of the barred galaxies with intrinsic ellipticities greater than 0.5 



Bar in disk galaxies 17 




Figure Al. The left panel plots the distribution of input (red histogram) and output (black, hatched histogram) values. In the right panel, the black 

line shows the correlation between output conce ntration and the bulge-to-disk ratio. The blue lines shows the relation between concentration and bulge-to-disk 
ratio from the 2-dimensional decompositions of Gadotti (2()kJ). 
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Figure A2. In the left panel, (ei ?arm „ pul — etarimrinsic) is plotted as a function of Rgo/Rso- The line shows the average value, with the error bars calculated using 
bootstrapping. In the middle panel, et armnpul (black line) and eb ar ,obs (red line) from real galaxies are plotted as a funciton of R90/R50. In the right panel, the 
distribution of simulated galaxies is plotted in the 2-dimensional plane of output concentration R90/R50 versus input ebar.intrinsic ■ The coloured contours denote 
the bar verification rate (F(bar)) as a function of position in the plane. The grey curve shows the average input ej,ar,mtrinsic that corresponds to an ei, arfiutout 
value of 0.5, as a function of Rgo/R^Q. 



